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O , Abstract. High Q^ NC and CC cross sections as measured at HERA can give 

information on two distinct areas of current interest. Firstly, supposing that all 
the electroweak parameters are well known, these cross-sections may be used to give 
K^ , information on parton distribution functions (PDFs) in the new kinematic regime at 

^ ' high X and high Q^. Secondly, supposing that PDFs are well known after evolution in 

Q^ from the kinematic regime where they are already measured, these cross-sections 
may be used to give information on electroweak parameters in a process where the 
exchanged boson is 'spacelike' rather than 'timelike'. WGl addressed itself to clarifying 
the limits of our present and possible future knowledge on both these points. 



1. Introduction 

The differential (Born) cross-section for charged lepton-nucleon scattering, mediated by 
the neutral current at high Q^ and accounting for the possible polarization of the lepton 
beam, is given by 

dxdq^ ~ Q^x 
where F2 is expressed in terms of parton distributions as 

F!f{x, Q2) = S,Af '^(g2) X (xg,(x, Q2) + ^q^^^,^ Q^))^ (2) 



r+Fr(a;,Q^)Tl^-xFf(x,Q^) , (1) 



and the sum is taken over flavours of quarks which are above threshold. (A contribution 
from the longitudinal structure function, F^, is neglected in Eq. |T]). Ai gives the 
couplings of the quarks and leptons to the currents 

A'^^''{Q^) = e] + 2e,ei{vi ± ai)v,Pz + {vi ± aif{v} + a1)Pl (3) 

These couplings depend on whether the polarization of the lepton beam is left (L, upper 
signs) or right (i?, lower signs) handed. The notation e; specifles the incoming lepton's 
charge such that e/ = —1. The vector and axial- vector couplings of the fermions are 
given by 

Vf = (T-if - 2e/ sin^ Ow), af = Tsf (4) 

where the deflnition holds good for any fermion, whether lepton or quark; T^f is the weak 
isospin, and 6w is the Weinberg angle 0. The term Pz accounts for the Z propagator 

Thus one can identify the contributions of 7 exhange, 7 — Z interference and Z exchange 
in the coupling Ai. Finally the parity violating structure function xF^ is given by 

xFf (x, g2) = S,i?/"^(Q2) X {xq,{x, Q2) _ xq^^x, Q^)), (6) 

where 

fif '^(g2) = ^2eiei{vi ± a{)a,Pz ± 2{vi ± aifv^a^Pl (7) 

The corresponding cross-sections for antilepton scattering are given by swapping L —^ R, 
i? — > L in the expressions for F2 and xF^ given in Eq. |^ and Eq. ^, and by substituting 
the antilepton charge e; = +1. 

The differential (Born) cross-sections for charged lepton-nucleon scattering, 
mediated by the charged current are given by 



d'a^^il^N) _ Gl M^ 



Y^F2ix,Q')TY-xFs{x,Q' 



dxdQ^ 4:7ix (g2 + M2,)2 

Expressing the structure functions in terms of parton distributions and accounting for 
polarization of the lepton beam gives 

whereas for antilepton scattering we have 

where the sums contain only the appropriate quarks or antiquarks for the charge of the 
current and the polarization of the lepton beam, P = {Nr — NL)/{Nfi + A''^,). 

t Neutrinos and charged leptons of the same family form weak isospin doublets with T3 = 1/2, —1/2 
respectively; and the quarks form similar weak isospin doublets, within the families (m, d), (c, s), (t, b), 
with T3 — 1/2, —1/2 respectively. For antiparticles both T3/ and e/ change sign. 



Clearly measurements of these cross-sections can give information both on 
electroweak parameters and on parton distribution functions (PDFs) . In Sec.^ and Sec.^ 
we consider how accurately we really know PDFs based on lower Q^ measurements, and 
hence how much HERA high Q^ measurements will improve our knowledge. Searches 
for new physics rely crucially on accurate determinations of event rates or cross-sections 
from standard model (SM) processes. In lepton-hadron or hadron-hadron collisions, this 
translates into a need to estimate uncertainties on the PDFs of the incident hadron. An 
example of such a need was seen when in 1997 both HERA experiments observed a slight 
excess of events at very high Q^ [|l] . Uncertainties on the SM cross-sections have been 
obtained by varying the input PDFs in the Monte-Carlo models. This method may not 
reflect the true PDF uncertainty, but rather highlight the different approaches used by 
global fitters in their PDF extractions. More reliable estimates of the uncertainties on 
PDFs are considered in Sec. 0. 

Recent work on parton distributions functions (PDF's) in the nucleon has focussed 
on probing the sea and gluon distribution at small x. The valence quark distributions 
have been considered to be relatively well understood. However, this is not really true 
at X >, 0.5. In particular our knowledge of the d quark distribution at high x is extracted 
from data taken on deuteron targets, and nuclear binding effects in the deuteron have 
not always been properly accounted. Our knowledge of the u quark density at high x 
has also been questioned by Kuhlmann et al, who recently proposed a toy model [^ 
which included the possibility of an additional contribution to the u quark distribution 
(beyond x > 0.75) as an explanation for both the initial HERA high Q^ anomaly [Q] 
and for the jet excess at high-P^ at CDF Q]. A precise knowledge of the u and d 
quark distribution at high x is very important at collider energies in order to estimate 
backgrounds to signals for new physics accurately. In addition, the value of d/u as a; — > 1 
is of theoretical interest in its own right. It is clear from Eqs. ^ and 10 that HERA CC 
data can greatly improve our knowledge of these high x valence quark densities. We 
discuss the current status of our knowledge in Sec. Ij. 

In Sec. H and Sec. || we discuss the possibilty of using high Q^ HERA data to gain 
information on electroweak parameters. We have specified the Born cross-sections for 
the processes above but if we aim for an eventual precision of 1 % then electroweak 
radiative corrections of order 0(a) have to be taken into account. We discuss the 
current status of implementation of the full electroweak radiative corrections and the 
implications for fitting electroweak parameters in Sec. |[ 

Measurements from LEP and the Tevatron tightly constrain the electroweak sector 
of the Standard Model (SM). In particular, the mass of the VT-boson, Mw, is known 
precisely from measurements of the properties of real, 'time- like', VT-bosons. In the 
SM description of charged current (CC) deep inelastic scattering (DIS) this same mass 
appears in the 'propagator' of the exchanged 'space-like' H^-boson (see Eq. ||). Hence, for 
the SM to be self-consistent the same value of Mw must be obtained from measurements 
of the space-like 'propagator mass' in CC DIS as is obtained from LEP and the Tevatron; 
any difference indicates new physics. In Sec. |] we discuss the present and future 



achievable accuracy of fits to M^ and G^ using HERA data. 

The workshop necessarily concerned itself with results which can be achieved at 
present or in the near future. In Sec. ^ we give a summary and discuss prospects for 
the more distant future. 

2. Uncertainties on Parton Distributions 

2.1. Experimental Uncertainties 

Measurements of cross-sections and structure functions inevitably involve detailed 
estimates of systematic uncertainties which are usually published with the 
measurements. The systematic uncertainties may be classified into two types: 
uncorrelated systematics in which the uncertainties fluctuate point-to-point (for example 
limitations on the statistics generated in Monte Carlo simulations of systematic effects) 
and correlated systematics which affect the data points in a coherent manner (for 
example, uncertainty in the luminosity determination). These systematic uncertainties 
may be propagated into the QCD fits used to extract parton densities and thereby allow 
an estimate to be made of experimental uncertainties in the extracted PDFs. A detailed 
derivation of the formalism necessary to make such estimates is given in |Q and both 
the HI and ZEUS collaborations have propagated systematic errors into their estimates 
of the errors on the gluon density at low a; ^, ^ . These experimental uncertainties are 
not the focus of the present section. They are thoroughly discussed by Botje 0. 

2.2. 'Theoretical' Uncertainties 

In addition to experimental uncertainties there are uncertainties which are of a more 
theoretical nature. A list of such sources of uncertainty is given below. 

• cXs- The size of as governs the NLO Q^ evolution of the PDFs. The Particle Data 
Group estimates an uncertainty of ±0.002 at a scale of M| [|]. 

• DGLAP solution methods. The DGLAP equations may be solved numerically 
using two classes of methods, namely x— space methods and n— space methods. 
Differences between these two approaches are generally small, but can be as large 
as 7% for the gluon density at low x @ 

• Heavy Quark Treatment. The way that heavy quarks are treated in DGLAP 
evolution can lead to uncertainties. In one approach the heavy quarks (c,6) are 
treated as massless partons and are evolved only above Q^ thresholds which are 
usually set equal to the masses squared. This method is not expected to give 
an accurate description of the data in the region of low Q^, where threshold effects 
appear as In Q'^/M^. In this region heavy quarks are considered as generated by the 



boson-gluon fusion process and the method of |]T0| is used. This method treats the 
mass logarithms correctly, but fails to match the asymptotic limit when Q^ ^ M^. 



The recent work of Lai et. al |Tl|, Martin et. al |]T2| and Buza et. al. [|T^ aims 



to have a smooth transition from the threshold region to the asymptotic hmit, but 
such an approach has not yet been implemented in the fitting programs used by 
the experimental collaborations. 

Target mass corrections. The identification of x with the fraction of the 
nucleon's momentum taken by the struck quark cannot be maintained when 
Q^ ~ M ^ (where M is the nucleon mass) and corrections to the formulae of the 
form x^ ■ M'^/Q^ are necessary 0]. 

Higher Twist effects. More complicated interactions involving more than one 
parton, like secondary interactions of quarks with the proton remnant, give rise to 
additional contributions to the structure functions. These higher twist terms are 
classified according to their Q^ dependence l/Q^", where n = is leading twist 
(twist=2), n=l is twist=4, etc. They are usually only important at high x [^ 
although higher twist effects at very low x have been considered recently [^ 



• Nuclear binding effects. Fixed target DIS experiments obtain information on 
the d density by scattering off a deuteron target. Theoretically the deuteron is 
often treated as a free neutron and proton; binding effects are neglected. Recently 
much work has been done in determining more realistic models of the deuteron. 
For more details see Sec. |[ 

Estimates of these uncertainties may be obtained by repeating several fits varying 
the conditions of each fit in order to ascertain the influence of the assumptions. In the 
following such a procedure is described with the aim of estimating uncertainties in the 
region of high x and high Q^, relevant to the new HERA data. 

2.3. Next-to-Leading Order QCD Fit of NC Data 



The DGLAP evolution equations [T^ are solved in the next-to-leading order (NLO) MS 
factorisation scheme using the program QCDNUM |18[. To ensure that the perturbative 



approximation is valid a starting scale of Qg = 4 GeV^ is taken at which four parton 
densities are parameterized. These are the up and down valence quarks [xu^ and 
xdy), the gluon (xg), and the sea quark distribution (xS). The sea quark density 
is assumed to have equal components from quarks and anti-quarks and is defined as 
xS = 2x(ia + d + s + c) at Qo where s = u/2 and d = u. Since the kinematic region of 
interest is at high Q^, a massless approach for heavy quarks is used. The xb density is 
evolved from zero at a Q^ threshold defined such that xb{x, Q"^) = for Q^< (5)^ GeV^ 
and the charm component is normalised to 2% of the sea quark density at Q^ = Qq, 
which gives a good description of the HI measurements [0 of ^2^^^, the charm induced 
structure function. 

The functional forms of the parton densities at Q^ are as follows: 

xu^{x,Ql) = A„x^"(l - x)'^"(l + D^x + E^y/x) 
xd^{x, Ql) = AdX^^{l - xf^il + DdX + EdVx) 
xg{x, Ql) = AgX^'{l - x)^«(l + DgX + Eg^) 



xS(x, Ql) = Asx^'^il - xf^{l + Dsx + Es^) 

The parameters A^, and Aj_ are determined by applying the valence counting 
rules and Ag is determined in terms of the other normalization parameters through 
the momentum sum rule. The strong coupling constant, a^, was fixed at the value 
«,(M|) = 0.118. 

Measurements of F2 from NMC [^ and BCDMS [^ on both proton and deuteron 
targets were used in the fit. In addition the F2 data from HI p were included for 
Q^ < 120 GeV^. Since there is no data to constrain the low x behaviour of valence 
quarks the parameters B^ and B^ were fixed in the fits to the values given in the 
MRSR2 parameterization. The infiuence of higher twist effects and deuteron binding 
effects was reduced by requiring the data to have Q^ > 10 GeV^, W'^ > 20 GeV^ as 
well as Q^ > 20 GeV^ for x > 0.5. The fixed target data were corrected for target 
mass corrections using the Georgi-Politzer approach [|14|, and deuteron binding effects 
as described in [^. Data with x > 0.7 were discarded, since in this region target 



mass corrections and deuteron binding effects are considered to be too large for reliable 
correction. The normalizations of all data sets were left free within the quoted luminosity 
uncertainty. The fits were performed using the program MINUIT [^ which minimised 



the x^ defined with statistical and uncorrected systematic errors added in quadrature. 
Several fits were performed varying in turn the input assumptions on the fit. 

• Og was varied by ±0.003 at the scale M|. 

• The infiuence of remaining higher twist effects was determined by increasing the 
Q2 cut from 10 GeV^ to 15 GeV^. 

• An uncertainty on the charm density was obtained by varying the the charm 
momentum fraction by factors of 2. 

• The strange component of the sea was varied by ±25%. 

• No deuteron corrections were applied to the data. 

The results of the fit are shown in the form of the reduced NC cross section (this 
term is defined in Ref. [^) in Fig |l|. Note that the uncertainty comes purely from 
the variation in the fitted PDFs, all electroweak parameters have been set to standard 
values. The total uncertainty on the NC cross-section is at the level of 1% increasing 
to 5% at a; = 0.65. This increase is largely due to the a^ variation and the Q^in cut 
applied to the data. The variation of the quark sea components has very little effect in 
this high Q^ region, however, the effect of the deuteron corrections can be as large as 
1%. The present study has considered only 'theoretical' uncertainties. It is interesting 
to compare our results with the analysis of Botje 0, where experimental systematics 
are considered (as well as an a^ variation). Botje finds the uncertainty on ai\ic to be 
4% rising to 11% at x = 0.4 and 5% to 8% at x = 0.6. Thus it seems that at lower 
Q^ the uncertainties are dominated by experimental errors, whereas at higher Q^ and x 
theoretical and experimental uncertainties are of similar size. 
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Figure 1. (Tjvc(e"^) extrapolated to high Q^. The band represents the uncertainties 
arising from variations in the input assumptions of the fit to NMC, BCDMS and HI 
data. 



3. Parton Distributions, d/u, and HERA CC Data 

3.1. Extraction of d/u at high x 

The most accurate information about valence quarks originates from fixed target 
electro/muoproduction data on proton and deuteron targets. The u valence quark 
distribution at high x is relatively well constrained by the proton structure function 
Ff. However, the d valence quark at high x is constrained by the neutron structure 
function F2, which is actually extracted from deuteron data. Therefore, there is an 
uncertainty in the d valence quark distribution from the corrections for nuclear binding 
effects in the deuteron. In past extractions of F^ from deuteron data, only Fermi motion 
corrections were considered, and other binding effects were assumed to be negligible. 
Recently, corrections for nuclear binding effects in the deuteron, F^jF^^^ ., have been 
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Figure 2. [a] The total correction for nuclear effects (binding and Fermi motion) in the 
deuteron, F2 / F^ , as a function of x, extracted from fits to the nuclear dependence 
of SLAC F2 electron scattering data, [b] Comparison of NMC F^^^ / F^ (corrected for 
nuclear effects) and the prediction in NLO using the MRS(R2) PDF with and without 
the proposed modification to the d/u ratio. From ref. ||2J| . 



extracted empirically from fits to the nuclear dependence of electron scattering data 
from SLAC experiments E139/140 [^. The empirical extraction uses a model proposed 
by Frankfurt and Strikman |^2|] in which all binding effects in the deuteron and heavy 
nuclear targets are assumed to scale with the nuclear density. The correction extracted in 
this empirical way is also in agreement (for x < 0.75) with recent theoretical calculations 
| 27| of nuclear binding effects in the deuteron. The total correction for nuclear binding 
in the deuteron is about 4% at x = 0.7 (shown in Fig. |^[a]), and in a direction which 
is opposite to that expected from the previous models which only included the Fermi 
motion effects. 

The ratio -F^/-Ff is directly related to d/u. In leading order QCD, 2F^/Ff — 1 ~ 
(1 + 4(i/«)/(4 + d/u) at high x. Bodek and Yang |2^ have recently performed an 
NLO analysis on the precise NMC F2/F2 data [^ taking deuteron binding corrections 
into account. The ratio ¥2^"^ /F2 is extracted by applying the nuclear binding correction 
F2 / ¥2^°^ to the F2/F2 data. As shown in Fig. ||[b], the standard PDF's do not describe 
this corrected Ff"*""//^!'. Since the u distribution is relatively well constrained, Bodek 
and Yang suggest a correction term to d/u in the standard PDF's (as a function of x), 
which changes the d distribution to fit the data. This correction term is parametrized as 
a simple quadratic form, 5{d/u) = (0.1±0.01)(x+ l)x for the MRS(R2) PDF, where the 
corrected d/u ratio is {d/u)' = {d/u) + 6{d/u). Based on this correction, an MRS(R2)- 
modified PDF is obtained, as shown in Fig 0[b]. The correction to other PDF's such 
as CTEQ3M is similar. Fig ^[a] illustrates the effect of this modification on the d/u 
ratio, where we can clearly see that d/u in the standard PDF's approaches as a; — > 1, 
whereas the modified d/u ratio approaches 0.2 ± 0.02 as x ^ 1, in agreement with a 
QCD prediction ||2^. Information on d/u which is free from nuclear effects may also be 
extracted from up/Vp data, but unfortunately this is rather inaccurate. Fig. ^[a] shows 
that the CDHSW |3^ data favour the modified PDF's at high x. Further information 
on d/u, which is not affected by the corrections for nuclear effects in the deuteron. 
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Figure 3. [a] The d/u distributions at Q^=16 GeV^ as a function of x for the standard 
and modified MRS(R2) PDF compared to the CDHSW data, [b] Comparison of the 
CDF W asymmetry data with NLO standard CTEQ3M, MRS(R2), and modified 
MRS(R2) as a function of the lepton rapidity. The standard CTEQ3M with a 
resummation calculation is also shown for comparison. From ref [ P5[ . 
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Figure 4. [a] The HERA charged current cross section data and [b] the CDF and DO 
inclusive jet cross section data are compared with both standard and modified PDF's. 
From ref I2I. 



can be extracted from W production data in hadron colliders. Fig ^[b] shows that the 
predicted W asymmetry calculated with the DYRAD NLO QCD program using the 
modified PDF is in much better agreement with recent CDF data [^] at large rapidity 
than standard PDF's. 

Clearly future e~^ /e~ HERA CC data can provide new information on the d/u ratio, 
free from nuclear effects. When the modified PDF at Q'^=16 GeV^ is evolved to Q^=10'^ 
GeV^ using the DGLAP NLO equations, it is found that the modified d distribution 
at X = 0.5 is increased by about 40 % in comparison to the standard d distribution. 
Fig. ^[a] shows the HERA e"*" CC cross section data compared to the predictions of 
the CTEQ4D distributions both with and without the modification. One can easily see 
that an increased luminosity of such data will be crucial in deciding the need for any 
modification. Fig. ^[b] shows that the modified PDF's also lead to an increase of 10% in 



the QCD prediction for the production rate of very high Pt jets p3[ in hadron colliders 
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Figure 5. Comparison of the modified NLO MRS(R2) PDF plus target mass and 
renormalon higher twist corrections, with low Q^ data, [a] Comparison of F2 and NLO 
prediction with and without higher twist contributions, [b] Comparison of R and NLO 
prediction with and without the higher twist contributions. From ref. pq]. 



3.2. Target mass, higher twist and PDFs at very high x 

Since all the standard PDF's, including the modified versions, are fit to data with x 
less than 0.75, Bodek and Yang also investigate the validity of the modified MRS(R2) 
at very high x (0.75 < x < 0.98) by comparing to Ff data at SLAC. Although the 
SLAC data at very high x are at reasonable values of Q^ {7 < Q^ < 31 GeV^), they are 
in a region in which non-perturbative effects such as target mass and higher twist are 
very large. The Georgi-Politzer calculation is used IQ for the target mass corrections 
(TM). These involve using the scaling variable 1^ 



2a;/(l + ^1 + 4M2xVQ2) instead of 
X. Since a complete calculation of higher twist effects is not available, data at lower Q^ 
{1 < Q^ < 7 GeV^ and x < 0.75) are used to obtain information on the size of these 



terms. Two approaches are used: an empirical method and the renormalon model ||34 
full details of the formalisms are given in 



5[. Both approaches describe the data well. 
Fig. 1^ illustrates the agreement for the renormalon higher twist approach. Thus one has 
confidence in extrapolating these approaches to investigate the very high x region. 

There is a wealth of SLAC data [^ in the region up to x = 0.98 and intermediate 
Q^ (7 < Q^ < 31 GeV^). Previous PDF fits have not used these data. The data for 
0.9 > X > 0.75 is in the DIS region, and the data for x > 0.9 is in the resonance region. 
It is worthwhile investigating the resonance region because duality arguments [^ 
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Figure 6. Comparison of SLAC F2 data witli the predictions of the modified 
MRS(R2), CTEQ4M and the CTEQ toy model at high x and 20 < Q^ < 31 GeV^. 
[a] Ratio to pQCD, [b] ratio to pQCD with TM effects, and [c] ratio to pQCD with 
TM and higher twist effects. From ref. p5||. 



indicate that the average behaviour of the resonances and elastic peak should follow the 
DIS scaling limit curve. Fig. ^ shows the ratio of the SLAC data to the predictions of the 
modified MRS(R2) at Q"^ values (21 < Q^ < 30 GeV^) where the elastic contribution is 
negligible. Fig ^[c] shows that, with the inclusion of target mass and renormalon higher 
twist effects, the very high x data from SLAC is remarkably well described by the 
modified MRS(R2) right up to x = 0.98. Fig|^[c] also illustrates that a good description 
is achieved using the empirical estimate of higher twist effects. Fig. ^ also shows that the 
CTEQ Toy model (with an additional 0.5% component of u quarks beyond x > 0.75) 
overestimates the SLAC data by a factor of three at x = 0.9 (DIS region). From these 
comparisons, it is found that the SLAC F2 data do not support the CTEQ Toy model 
which proposed an additional u quark contribution at high x as an explanation of the 
initial HERA high Q^ anomaly and the CDF high-Pj jet excess. 

Thus it seems that there is a need for modification of current standard PDFs at high 
X, but this conclusion is dependent on the correctness of deuteron binding calculations. 
Future precise charged current data from HERA will provide vital information on valence 
quark densities at high x and high Q^, free from this uncertainty. 

4. Electroweak Radiative Corrections 



Electroweak radiative corrections of order 0{a) have to be taken into account if a 
precision of 1 % is to be achieved in physics analyses at HERA. Pure QED corrections 



contain contributions which are enhanced by logarithmic factors hke \og{Q^/ml). In 
addition, radiation of photons can shift kinematic variables from very large to very 
small values inducing additional enhancement factors for NC scattering. Furthermore, 
purely weak one-loop corrections will eventually be important for a measurement of M^ 
with a precision of 0(100 MeV). 

As discussed in much detail in the literature |^ , the differential cross section for NC 
scattering when measured in terms of electron variables receives QED corrections which 
are large and negative at small ye and large Xe, whereas at large i/e and small Xe large 
positive corrections may reach the level of 100 %. These corrections can be reduced, 
sometimes considerably, with measurements based on hadronic or mixed kinematic 
variables. Self energy corrections (the running of a(Q^)) are of order 10 % and, finally, 
purely weak corrections are negligible at low Q^, but increase to the level of a few per 
cent at very large Q^. 

In the following we give a more detailed discussion of electroweak radiative 



corrections to the charged current process [^, |3^, ^, firstly since their treatment 
is more involved and has been less well described in the literature than the NC case; 
secondly because the implementations in available numerical programs appear to have 
larger uncertainties than the corresponding NC codes. 

Fig. 1^ shows the complete set of electroweak one-loop Feynman diagrams 
contributing to eg/ -^ z/g// scattering. The list comprises corrections to the euW and 
QfQfiW vertices, self energy corrections to the external lines, graphs for the W self energy 
and box diagrams. Diagrams with an additional photon in the loop gives rise to infrared 
divergent contributions. They can be regularized, for example by using a finite mass for 
the photon, and cancel against similarly infrared divergent contributions from radiative 
CC scattering eg/ -^ ^Qfl (see Fig. |^) when calculating the cross section. 

The calculation is more complicated than the one for NC scattering since for the 
CC case there is no simple enough subset of Feynman diagrams which is gauge invariant. 
In particular, all four diagrams of Fig. |] describing the emission of a photon from the 
incoming lepton (a), from the incoming (b) and outgoing quark (c) and from the internal 
W (d) are necessary to obtain a physically meaningful result. 

Internal fermion propagators in the diagrams of Fig. |^ give rise to terms with pole 
factors which become large when the photon is collinear with one of the external charged 
fermions. The differential cross section for radiative CC scattering can thus be written 
as 



rfV 



= -— + h - — - + non — pole terms, (11) 

, kl kq kq' 



dxdQ'^d?k 

where the notation shown in Fig. ^(a) for the particle's 4-momenta was used. This 
shows that the numerical calculation of the 7-inclusive cross section has to cope with 
nearly-divergent behavior in different phase space regions. Note that the separation into 
pole factors in Eq. ^ is not unique: finite (i.e. non-pole) terms can be shifted freely 
between the three contributions. A specific separation of the radiative cross section 
like that given in Eq. |ll] is needed for the construction of a well-behaved numerical 




Figure 7. One-loop Feynman diagrams for eq -^ uq' scattering. 







(a) (b) (c) (d) 

Figure 8. Feynman diagrams for radiative charged current scattering eq — > vq'^ and 
notation of 4- momenta for external particles. 



algorithm, but the individual contributions do not have a well-defined physical meaning 
when taken separately. 

To obtain a gauge invariant, i.e. physically meaningful separation, one can organize 
the terms contributing to the complete cross section according to their dependence on 
the electric charge of the incoming particles. Replacing the charge of the W^ by that 
of the incoming lepton, e^ = cvk, and that of the final quark qji by that of the incoming 
quark g/ and the lepton, e// = e/ + e/, one can write 



dxdQ'^(Pk 



= 6; /lep + eiCfli^t + efl^^^. (12) 

The individual parts in this separation are gauge invariant and due to their association to 
fermion charges it is justified to call them "leptonic", "interference", and "quarkonic", as 
indicated by the subscripts in Eq. |1^. A closer inspection of the diagrams in Fig. ^ shows 
that Jicp contains pole terms oc 1/kl and 1/fcg', whereas /qua contains terms oc 1/kq and 
1/kq' . The interference contribution lint receives contributions from all three pole factors 



1/kl, 1/kq and 1/kq' . After integration over the photon momentum, the pole factors give 
rise to logarithms of the corresponding fermion masses: log(Q^/mg), log(Q^/'m^ ) and 
\og{Q'^ /m?q J. Here we assume that collinear singularities are regularized by introducing 
finite masses for quarks. 

The virtual corrections from Fig. ^ can be separated in the same way. After 
combining real and virtual corrections, not only do infrared divergent contributions 
cancel (allowing the removal of the regularization by the finite photon mass), but also 
many mass logarithms disappear in the complete final result: writing the complete 
0(a)-corrected cross section similarly to Eq. |T2l : 



dxdQ'^(Pk 



efJlep + eiCfJint + ej Jqua, (13) 



one finds that i) all logarithms related to the final quark disappear, i.e. Jiep contains 
log{Q^/ml) terms only, Jqua contains log(Q^/m^ .) terms only and ii) all mass logarithms 
cancel in Jjnt, i.e. the interference contribution contains no mass singularities. Therefore 
the effect of photonic corrections on the cross section for CC scattering can be 
summarized as follows: 

• Leptonic contributions proportional to ef contain terms which are enhanced by 
the large logarithm \og{Q'^/ml). This contribution has a dependence on kinematic 
variables which is very similar to the NC case except that the enhancement at large 
He is much less pronounced because of the much weaker dependence of the CC Born 
cross section on Q^. 

• Interference contributions are of order 0(a/7r) without logarithmic enhancement. 

• Quarkonic contributions contain uncancelled log{Q'^/m'^,) terms. These mass 
singularities have the same structure as those appearing in a calculation of 0{as) 
QCD corrections due to the emission of gluons from incoming and scattered quarks 
(in dimensional regularization one would find the same 1/e divergences apart from 
differing charge factors: ael instead of asCp). As with QCD corrections, the 
mass singularities due to QED corrections can be factorized and absorbed into 
renormalized parton distribution functions. Thereby the final results do not depend 
any more on a regularization parameter (like the quark mass), but the Q^-evolution 
equations for parton distribution functions receive an additional QED term: 

Q^J.,,^^, Q.) ^ ^ ^' f (C.P„,(.),, (f . 0^) + P,,(.), (f , Q 



Numerical estimates [Q show that such effects reach the level of 1 % only at very 
large x and Q^. 

One should note that the available numerical programs either simply ignore the 
quarkonic contributions or suppress them by using a relatively large value for the quark 
mass (of the order of several GeV). Another level of approximation justified by the 



above discussion is to neglect the interference contribution. Numerical results for the 
individual contributions have been given in |3^, ^ . 



Having separated the purely photonic corrections as described above, there still 
remain purely weak contributions (most prominently, but not exclusively, due to the 
W self energy) which can be combined into over-all form factors [^, 0. For electron 
scattering with unpolarized beams Eq. H becomes: 



da 



dxdQ"^ 



e~p \^ ' ^''^W / K f=u,c 



+(i-2/)'(pcc(^,Q'))' E <iA^^Q')\- (15) 

f=d,s ) 

The form factors p^(^(a;, Q'^) depend on Q'^ and x as well as on the type of the scattered 
quark (due to the presence of box diagrams). The above expression assumes that the 
CC cross section is normalized with the help of the /i-decay constant G^. 

The cross section also depends explicitly on the W boson mass and, through the 
form factors Pcc^ o^ ^^^ other parameters of the electroweak standard model. 

4.1. Choice of independent parameters 

One possible set of independent parameters is a, nit (the top mass), Mh (the Higgs 
mass), Mw, Mz (in addition, there are light fermion masses and CKM matrix elements, 
which we do not consider). A special role is played by a which is kept fixed since it 
is well-measured and is a parameter of QED which is embedded in the full electroweak 
theory. The Higgs mass is also special: because of the weak dependence on Mh in 
present day experiments, one can keep it as an external parameter and investigate the 
effect of changing it separately (see Sec. |). Thus the essential parameters are Mw, Mz 
and rrit, since G^ is related to these parameters by the SM expression given in Eq. 16. 
We consider the following two possible scenarios: 

• Choose Gfj,, Mw and Mz as independent parameters. Then ruf is a prediction and 
a specific choice of values for G^, Mw and Mz will predict a value for rrit (possibly 
in disagreement with experimental results). Fixing G^ and Mz to their measured 
values (from jS decay and LEP experiments, respectively) the CC cross section 
depends on Mw and a measurement at HERA can be interpreted as a measurement 
of the W mass. Note that this assumes the validity of the Standard Model and 
the measurement is therefore not a direct PF-mass measurement. (This is also 
true for scenario 2). Alternatively, keeping Mw and Mz fixed, the charged current 
cross section could be interpreted as a measurement of G^. This is interesting 
since it allows a direct comparison of the high-Q^ experiments at HERA with the 
low-energy measurement of p decay 

• Choose Mw, Mz and rrit as independent parameters (on-shell scheme). Then G^ 



is a prediction of the theory given by 



Tia 1 . , 2 . ^M 



g.= ^:'"„ '., with 4 = 1-T7f (16) 



and Ar, which embodies the radiative corrections to the /i-decay, is a function of 
Mw, Mz, rrit (and a, Mh, of course). Ar has to be taken into account in a consistent 
0{a) prescription and should be combined with the other 0{a) corrections in the 
CC cross section formula, i.e. replacing 

Because of the inclusion of Ar the new form factors have a strong dependence on 
rrit since 

Ar = Aa - ^Ap + Ar'"'='" with Ap = ^' * . (18) 

Furthermore the dependence of the CC cross section on M^y is much stronger due 
to the replacement 

( <■ V^f ^ ^V (19) 

Applied to the situation at HERA, keeping Mz and rrit at their experimental values, 
the measurement of the CC cross section can once more be interpreted as an indirect 



measurement of M]y p3| , |4^ which is, through Eq. |1^, equivalent to a measurement 
of Gf^, however in a process at large momentum transfer. 

Both schemes are theoretically equivalent, a specific choice should be motivated by the 
aim of the analysis. We repeat: the above discussion focussed on the essential parameters 
Mw, Mz, rrit, but, in addition, a weak dependence on Mh is always present, both in 
Ar as well as in the form factors p^C' ^^^^ ^^ quantified further in Sec. |^ 

It might be instructive to see how the situation simplifies in the case when radiative 
corrections are neglected altogether. Then rrit and Mh do not play any role for CC DIS. 
In scenario 1 where G^ and M^ are considered as independent parameters, Mz is also 
irrelevant. With fixed G^, the measured cross section can be fitted by varying Mw The 
resulting value has been called a 'propagator mass' in previous experimental literature 
since it is due to the dependence of the cross section on Mw via the propagator only. 
Keeping Mw fixed leaves no possibility to adjust the Q^ dependence, and then the 
measurement only determines the normalization (in terms of G^) of the cross section. 
Scenario 2, with fixed Mz, allows us to combine the sensitivity to the propagator and 
to the normalization in one common parameter, i.e. the W mass, as made explicit in 
Eq. |l^. In Sec. |^ we consider both scenarios. 

4-2. Numerical results and discussion of uncertainties 

The presently available numerical programs for the calculation of the CC cross section 
do not all take into account the complete set of one-loop and one-photon corrections 
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Figure 9. Comparison of radiative corrections for the charged current cross section, 
6cc = da/dxidyi\Q,^s / da/dxidyi\^^^^ — 1, as a function of yi for fixed values of Xi. 
The full curves include the complete QED and weak corrections, the dashed curves are 
without purely weak corrections and the dotted lines show the leptonic QED corrections 
only (all from epcctot^. Open symbols are obtained from DJANGOH and include leptonic 
QED and purely weak corrections. 



described above. Two programs with the complete set of 0{a) corrections (DISEPW 
38| and epcctot ^^) have been compared in the HERA workshop [^ and found to 



agree well except in the phase space corners of small x, small y or large x, all y. These 
programs are, however, not very well-suited for an application to realistic experiments 
since they do not allow to apply experimental cuts and, secondly, since they are restricted 
to the use of leptonic variables whereas experiments have to use hadronic variables. The 
Monte Carlo event generator HERACLES/DJANGOH [Q circumvents these two restrictions, 
however, it has the CC radiative corrections implemented in an approximation where 
photonic radiative corrections from lepton-quark interference and from quark radiation 
are neglected. Theoretically this is motivated to the extent that an accuracy of 1% is 
sufficient, as discussed above. 

Fig. 1^ shows a comparison of results obtained with the most recent version of 
DJANGOH with the complete 0{a) results from epcctot for lepton variables. At small 
values of Xe, the comparison is satisfactory, but at larger values of Xe, the Monte Carlo 
results lie above those from epcctot. Here improvements are still needed. More 
detailed comparisons, in particular between DJANGOH and HECTOR |^, are presently 
being performed [^]. 

Apart from the approximations in available program implementations, there are 
additional intrinsic uncertainties of predictions for the CC cross section. First, a 



theoretically consistent treatment of quarkonic radiation as described above is not 
available yet. The uncertainty from this source is most likely below 1%, and is thus 
negligible except at very large Q^ |^. Secondly, higher-order corrections are to be 
expected, first of all from multi-photon effects, but also from combined QED-QCD 
corrections of O (acts) and possibly from two- loop weak corrections. Multi-photon effects 
due to soft photons can be treated by exponentiation. This effect can be taken into 
account only in epcctot or HECTOR. It is expected to be numerically significant at low 
?/e, low Xe (0(2%) for Xe < 0.1) and very large Xe (0(10%) for Xe > 0.9) Q, i.e. 
close to the phase space boundaries. Two-photon contributions of O(a^) in the leading 
logarithmic approximation are also known and can be calculated with the help of HECTOR 
but no numerical evaluation of this contribution has been performed for CC scattering. 
0{aas) corrections are not expected to become important. An indication of this comes 
from the fact that in NC scattering the larger 0{as) correction to F3 stays below 1.5%, 
for values of x and Q^ relevant at HERA. Two-loop weak contributions are also not 
expected to be important for the cross section itself. However, two-loop effects can 
affect the value of M\y at the level of several tens of MeV when determined from G^ 
via Eq. |16|. The dominating higher-than one-loop effects to the G^ — My/ relation are 
implemented in all availabe programs, but not all of the many small effects discussed 
above have been taken into account while the dominating source of uncertainty is still 
the statistical error. 

5. The Determination of M^/ in Charged Current Deep Inelatic Scattering 

5.1. Introduction 

Recently the ZEUS and HI experiments at the HERA positron-proton collider 
announced preliminary measurements of the single and double differential cross sections 



for the CC reaction e'^P -^ u^X |^, Q. The measurements were made using data 
sets corresponding to a luminosity of ~ 40 pb^^ per experiment. The CC data 
is now sufficiently precise to make it timely to re-examine the extent to which the 
HERA experiments can contribute to constraints on the electroweak sector of the SM. 
In the following the statistical precision with which the parameters Mw or G^j may 
be determined will be investigated. Hence, it is the dependence of the CC cross 
section on G^ or M\y which is of primary importance and small differences between 
different approaches to the calculation of the radiative corrections themselves will not 
be discussed. There is no attempt in what follows to discuss the experimental systematic 
uncertainties involved in measurements such as those proposed. 

5.2. Sensitivity of CC DIS to Mw 

The Born cross-section for the reaction e^P -^ UgX has been given in Eq. ^. In this 
expression the coupling of the W^-boson to the fermions is contained in G^ while the 
VT-boson propagator contributes the factor [M^r/{M^, + Q^)]"^- 



The mass of the exchanged boson can be determined from the Q"^ dependence of 
the CC cross section since the propagator factor ([M^/(M^ + Q^)]^) depends exphcitly 
on Mw- This is the 'propagator' mass mentioned in Sec. ^ (scenario 1 with fixed G^). 
Both the ZEUS and HI collaborations have performed such a determination with the 
results g|, g 

Mw = 78.Qtli [stat] tH [syst] t\i [PDF] GeV (ZEUS), (20) 

Mw = 81.2tli [stat]+4:3 [syst] +|^ [PDF] GeV (HI). (21) 

The uncertainty on the value of Mw arising from the PDFs is quoted separately above. 
This error is determined by performing the fit using a variety of PDFs (additionally, the 
error quoted by HI includes a 4% error on the CC radiative corrections). 

The data may also be used to determine G^ or Mw (assuming scenario 2). The 
precision with which G^ can be determined is directly related to the precision with which 
the cross section is measured. The size of the e~^P data set available to ZEUS and HI 
yields ~ 1000 CC events per experiment. Hence the statistical uncertainty on such a 
direct determination of G^ will be ~ 1.7%. Note that at HERA such a determination 
of G^ will be performed at an effective Q^ in excess of ~ 400 GeV^ and hence will be 
complementary to the precise measurement of G^ obtained from muon decay where the 
momentum scale is set by the mass of the muon. 

At O (a) in the SM G^ depends upon Mw through Eq. |16[ In the following we shall 



set a, Mz, and all fermion masses, other than rrit to the values quoted in reference |51 
As shown in Fig |10|, G^ depends strongly on Mw and less strongly on the mass of the 
Higgs boson, M^, or the mass of the top quark, rrit. Hence, within the context of the 
SM, the greatest sensitivity to Mw may be obtained by combining the Mw dependence 
of the propagator term in the CC cross section with the Mw dependence of G^. A naive 
propagation of the statistical error from a sample of 1000 CC events leads to an error 
of ~ 300 MeV on the parameter Mw- More detailed numerical estimates are given in 
the next section. The figure shows the expected statistical error on a measurement of 
Gfj^. It can be seen that the variation of G^ with nit or Mh at fixed Mw is negligible 
compared to the expected statistical error. 

Before proceeding to numerical estimates of the precision with which Mw or G^ 
can be extracted from present CC data it is necessary to consider the dependence of the 
radiative corrections to the Born cross section (Eq. §) on Mw- The radiative corrections 
to the Born cross section given in Eq. |l^ may be formally expressed in terms of a 
correction factor, ^Rad, defined by the equation 

" "^Rad _ " '-'Born n , r 1 ,'00^ 

The full O (a) radiative corrections may be split into two pieces by writing 

1 + ^Rad = {1 + 6^ + 5w) . (23) 

A real photon may be radiated from one of the particles involved in the hard scattering 
or appear as a virtual particle in a one-loop diagram. Such corrections contribute to the 
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Figure 10. The dependence of G^ on Mw embodied by the O {a) expression of 
equation |l^. The thick sohd line shows the resuh when nit is set to 175 GeV and Mh 
to 100 GeV. The thin sohd, dashed, dotted and dash dotted hnes show the dependence 
of G^ on Mw for the values of mt and Mh indicated. The point with error bar shows 
the expected statistical uncertainty on G^. 



term 5^. All other O (a) contributions are included in the term 5w. The size of these 



corrections is shown in Fig 11. The full radiative correction, ^Rad? was calculated using 



the program HECTOR [0, while 5w was calculated using the program EPRC |^. It 



can be seen from Fig [TT| that (5Rad can be as large as ±20% while ^w is never larger than 
-3%. 

The radiative corrections to CC DIS depend on G^; i.e. 6 = 6 (Gfj), where 6 may 
be either (5Rad or 5w- Thus for the determination of G^ (or M^i/) from the Born cross 
section to be valid it is important that (^Rad should depend only weakly on G^. The 
sensitivity of 6-Rad and 6^ to G^ is shown in Fig ^ where the ratio 



R 



6{g, 



1.226 X 10-5 GeV" 



6{G, 



1.166 X 10-5 GeV" 



(24) 



is plotted as a function of x for several fixed values of Q^. The reference value of G^ 
[Gfj_ = 1.166389 X 10~^ GeV~^) was taken from reference |^. Given a statistical error of 



~ 1.7% the value G^ = 1.226 x 10~^ GeV~^ corresponds to raising G^ by three standard 
deviations from the reference value. Fig |12| shows that the sensitivity of (5w to such a 
change in G^ is small (< 1%). The full radiative correction, (5Rad, shows a still weaker 
dependence on G^ since ^Rad is dominated by the photonic corrections contained in S^. 
Given the weak dependence of the radiative corrections on G^ it is reasonable to 
estimate the precision with which G^ (or M^y) may be extracted by making a fit to 
the Born cross section. In practice, once G^ has been determined from such a fit, 
the extraction of the experimental cross section would have to be iterated using the 
new value of G^ in order to verify that the results obtained were stable against such 
modification. 

5.3. Numerical Estimates 

The program HERACLES [^ was used to estimate the statistical error corresponding 
to a luminosity of 40 pb~^ of e^P data (a total of 1020 CC events) distributed among 
eight bins in Q^. The bins were equally spaced in log^g Q^ ^^ Q^ i^ the range 
400 < g2 < 40000 GeV2. Eq | was then used to evaluate da^gJdQ'^. The CTEQ4D 
PDFs were used through out. 

An estimate of the precision with which G^ may be determined was obtained as 
follows. The value of Mw was fixed at Mw = 80.41 GeV and a x^ fit of the differential 
cross section, da{G^)/dQ'^, performed. This is scenario 1 with M]y fixed, as described 
in Sec. §. The fit yielded an error on G^ of ±0.018 x 10~^ GeV^^ which corresponds 
to ±1.5%. A slightly improved error may be obtained by considering scenario 2 which 
combines the sensitivity to the propagator and to the normalization into one common 
parameter, using Eq. |16[ The error on G^ is then ±0.016 x 10^^ GeV^^ (±1.4%). This 
last result may be re-expressed to give an estimate of the error on the value of M^y, 
using Eq. |19|, as suggested in Sec. ^. Such a fit was performed and yielded a statistical 
error on the parameter M^ of ±0.24 GeV. These results agree well with similar fits 



reported in references [Q and Q . 

6. Summary and Outlook 

We now consider what we have learnt and what we can learn given various running 
conditions: e"^/e~ beams, polarization of the beam, beam energies. 

We have seen that HERA data has been used to measure the 'propagator' mass Mw 
and that the precision on such a measurement may be improved by using information 
on the normalization as well as the shape of the cross- sect ion. Future improvements 



in luminosity (~ lOOOpb^ ) may allow us to achieve an error of ±55MeV on M\y ||44 
Sensitivity is greatest with electrons rather than positrons for both NC and CC processes 
and NC data is most useful with left handed polarization of the lepton beam (which 
enhances the 7/^ terms, see Eqs. |l|- 0). Indeed 70% polarization is worth a factor of 4 
in luminosity. 
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Figure 11. Radiative corrections to charged current deep inelastic scattering. 
The solid line shows the full radiative correction, (Jpiad: calculated using the program 
HECTOR 1^. The dashed line shows the electroweak contribution to the radiative 
corrections, 6^^, calculated using the program EPRC Eaj. 
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Figure 12. The dependence of the radiative corrections to charged current deep 
inelastic scattering on the parameter G^. The ratio i?, defined in the text, shows the 
fractional change in the radiative correction as Gfj, is varied from 1.116x 10~5 GeV~^ to 
1.226 X 10~5 GeV^^. The solid line shows the sensitivity of the full radiative correction 
(evaluated using HECTOR [Q) to G^. The dashed line shows the sensitivity of the 
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One may also be able to measure the weak neutral current couplings of quarks, f„, 
ci-u, Vd, ttd (see Eqs. |I]- 1^, as explored in Ref. [0. With unpolarized e"*" and e~ beams 
ratios of NC and CC cross-sections can be used to determine Qu, whereas polarized 
beams allow the extraction of Vu- Measurements of NC cross-sections alone may be 
used in a fit to determine all four couplings. Polarization of the beams is essential 
to achieve reasonable precision. One needs ~ 250pb~^ luminosity for each of the four 
lepton beam charge/polarization combinations to make measurements with a precision 
of ~ 10%. 

We have seen that HERA CC data can improve our knowledge of the u and d valence 
distributions (Sec. |^). It seems that knowledge of the d quark is the most crucial since 
the u quark is already much better known, not being subject to uncertainties due to 
deuteron binding corrections. To improve our knowledge we will require further positron 
beam running. Since we would like to explore the valence distributions at the highest 
possible values of x it may also be worth lowering the proton beam energy to achieve 
higher x values at the same Q^ (x = Q'^/{ys)). This would have two further advantages: 
it would allow substantial overlap with fixed target data at lower Q^ and thus constrain 
normalizations, and it would allow a model independent measurement of the longitudinal 
structure function F^ ^5 . 



Can we extract other interesting structure functions or PDFs with future high 
luminosity, high Q^, NC and CC data? This has been explored in |Q. Polarization is 



not explicitly useful for this purpose so we will assume P = in the formulae (though 
they can be generalized to any value of P). If we express Eqs. ^ and ^ in terms of 
reduced cross-sections we obtain 

^cc{e~p) = [x{u + c) + (1 - yYx{d + s)] (25) 

^cc{e^p) = [xiu + c) + (1 - yfx{d + s)] . (26) 

Then we can define the sum and difference of these reduced cross-sections (x±{CC) = 
^cc{e~p) ± acc{e + p), so that 

a+{CC)=xU + {l-yYxD (27) 

a_{CC) = xu^ - (1 - y^xd^ (28) 

where U stands for the sum of the w-type quarks and antiquarks {u,u, c, c), D for the 
sum of the d-tjpe {d, d, s, s) and u^ and d^ stand for the valence quark densities. 
Similarly the sum and difference of the NC e~^p and e~p reduced cross-sections give 

a+{NC) = AuxU + AdxD = F2 (29) 

a_{NC) = B^xuy + Bdxd^ = xFs (30) 

where A^, -B„ are the same for all w-type quarks and antiquarks and Ad, Bd are the 
same for all (i-type quarks and antiquarks. 

Thus from the CC cross-sections we can mainly gain information on the valence 
distributions at high x, whereas from the NC cross-sections we may measure xF^ for 



the first time. We may also consider making combinations of all four cross-sections 
d^l_{NC /CC) in order to extract f/, D, u^, d^ for the complete x range. A detailed 
analysis of this possibility indicates that 200pb~^ per e^/e~ beam will be necessary to 
make measurements of ~ 10% accuracy [0]. Finally, can HERA shed light on the charm 
and strange components of the sea? Eqs. ^ ^ indicate how the CC processes may be 
studied at low x and high and low y to look for evidence of the appropriate sign of D 
or K hadrons in the final state. 

In summary it seems that running conditions with an equal partition between e"^/e— 
and left/right handed polarization will give the most fiexibilty for future physics, and 
that running at a reduced proton beam energy is of interest to physics at high Q^ as 
well as at low Q^. 
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